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Abstract 

Classical MD data on the charge-charge dynamic structure factor of two-component plasmas 
(TCP) modelled in Ref. jl] are analyzed using the sum rules and other exact relations. The 
convergent power moments of the imaginary part of the model system dielectric function are 
expressed in terms of its partial static structure factors, which are computed by the method of 
hypernetted-chains using the Deutsch effective potential. High-frequency asymptotic behavior of 
the dielectric function is specified to include the effects of inverse bremsstrahlung. The agreement 
with the MD data is improved, and important statistical characteristics of the model TCP, such 
as the probability to find both electron and ion at one point, are determined. 

PACS numbers: 52.27.Gr, 52.25.Mq, 52.27.Aj, 52.65.Yy 
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I. INTRODUCTION 



Classical molecular dynamics (MD) simulations of model two-component plasmas (TCP) 
were carried out in Ref. jl| more than twenty five years ago. The system modelled in this 
pioneer work was a fully ionized strongly coupled hydrogen TCP of temperature T consisting 
of electrons (a = e, = —1) with the number density n and protons (a = i, Zi = +1) with 
the same number density. The dynamic and static characteristics studied were the charge- 
charge dynamic structure factor, the species partial static structure factors and the static 
radial distribution functions, etc. Classical statistical averages were computed on the basis 
of the ergodic hypothesis while the quantum delocalization preventing the collapse were 
taken into account through the use of the Deutsch pair effective potential 2j arising from 
the quantum-diffraction effects without including the exchange or symmetry contribution. 



VJafe (?^) = ^a^b— (1 - exp (-Kafcr)) , a,b = e,i, (1) 



where Hab is the reduced mass of an a — 6 pair and [3~^ is the temperature in energy units. 
The fact that the potentials ([1]) remain finite as r — > is a consequence of the uncertainty 
principle and prevents the collapse to which we have already referred. In the temperature 
range of interest 

K.a»l, 

being the Wigner-Seitz or the "ion sphere" radius. Thus the effective ion-ion interaction is 
virtually identical to the bare Coulomb potential at all separations. 

In addition, in Ref. l| the potential ([T]) was also employed to determine the static 
properties in the hypernetted-chain (HNC) approximation. 

Since 1981, little effort was made to simulate this benchmark high-energy density system 
and study its static and dynamic properties. We mention a few related works: (i) the static 



electrical conductivity was studied in Refs. 



Sj and 4(], (ii) the TCP dynamic characteristics 



in a different range of values of the wavenumber were investigated in Ref. j^. 

The results of Ref. were analyzed in Ref. {g] using the sum rules and other exact 
relations. An overall agreement with the MD results was obtained in Ref. j^, where the 
frequency moments of the imaginary part of the plasma inverse dielectric function, e~^(/c, cu), 
(the sum rules) were calculated for the bare Coulomb potential. 



Our aim here is to reexamine the simulation data of Ref. ^ and the theoretical results 
of Ref. within the moment approach and using the method of effective potentials to 
evaluate the static characteristics of the TCP. 

Mathematical details of the moment approach are provided in the following Section. A 
model Nevanlinna parameter function taking into account the fractional asymptotic form of 
the imaginary part of the system dielectric function is suggested in Sect. IIIBj in Sect. Ill CI 
the convergent power moments of the loss function [— Im {k , uj) / u] of the model system 
are calculated using the Deutsch effective potential and are computed by the method of 
hypernetted-chains using the same potential; for a recent review of the method of effective 
potentials see Ref. Q]. Numerical results and conclusions are presented in Sects. Ill and IV, 
respectively. The agreement with the MD data is improved significantly, and, simultaneously, 
important statistical characteristics of the model TCP in concern, like the probability to find 
both electron and ion at one point [gei (0)], are determined, direct and exchange interaction 
contributions being compared. 

II. THE METHOD OF MOMENTS 
A. The background 

The MD results of Ref. jl] on the charge-charge dynamic structure factor Szz{k,uj) were 
modelled in Ref. [6] using the moment approach which automatically takes into account the 
sum rules and other exact relations. 

The starting point in the application of the method of moments to the calculation 
of the system dynamic correlation function is the fluctuation-dissipation theorem (FDT) 
which relates the latter to the system dissipation characteristic, the Green function, whose 
power frequency moments, by virtue of the Kubo theory of linear response, can be directly 
expressed in terms of the static correlators of the time derivatives of the system observables. 
These static characteristics can be expressed, using the system model Hamiltonian, in terms 
of the system structural static correlation functions like the radial distribution functions 
or the static structure factors, and this is the only step of our approach where the system 
model interferes. Otherwise the relations we use are model-free. They are based on mathe- 
matical results which are independent of the physical details, i.e., the interaction potential. 
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This permits to apply the moment approach to non-perturbative systems which lack small 
parameters, such as the one we deal with here. 

Indeed, in a classical Coulomb system the characteristic perturbation parameter is the 
potential energy of two charges at an average distance, i.e., the Wigner-Seitz radius a, 
divided by the characteristic kinetic energy or the system temperature: 

r = I3e^/a. 

The achievement of the work by Hansen and McDonald [l| was that the values of F > 1 
(precisely, F = 0.5 and F = 2) and = 0.4 and = 1, where 

a ame^ 



ao being the Bohr radius, were considered, which correspond to strongly coupled Coulomb 
systems or high-energy density (HED) matter. Under these conditions the Landau length, 
the Wigner-Seitz radius, and the Debye radius 

n 

A 



are of the same order of magnitude, so that no effective screening takes place: there is only 
one or less particles in the Debye sphere. Hence, standard theoretical treatments, e.g., the 
kinetic theory are not applicable, and we need alternative approaches such as the one based 
on the method of moments. 

There is an additional dimensionless parameter which is commonly used to characterize 
the kind of statistical systems we consider. This is the so-called degeneracy parameter, 

= (5Ef = 1.84159— 

which compares the Fermi energy of the system, Ep, to its thermal energy. For the ther- 
modynamic conditions considered in Ref. l| and here > 1, which means that strongly 
coupled two-component plasmas are intrinsically quantum statistical systems, and we might 
expect their physical properties to be greatly influenced by their quantum mechanical nature. 

Nevertheless, due to the significant difficulties encountered in the analytical and computa- 
tional modelling of these systems, it is a commonplace to use classical statistical approaches 
to investigate multicomponent plasmas. These classical approximations, like MD or HNC 
calculations, do require the application of potentials which effectively take into account the 



quantum mechanical nature of these systems, especially the quantum diffraction preventing 
the classical Coulomb collapse. 

As we have said, the only step where the physical model under consideration appears 
explicitly in the method of moments is the modelling of the system Hamiltonian required to 
compute the frequency moments of the Green function (see Sect. Ill C\i . Strictly speaking, 
these calculations can be carried out according to the rules of quantum statistical mechanics, 
and therefore no effective potential is needed. This means that the potential appearing in 
the Hamiltonian must be understood as the real interaction among charged particles at the 
microscopic level, which is, of course, the Coulomb potential. 

However, our main aim here is to compare the static and dynamic characteristics obtained 



for the TCP with those of Ref. [l]. Thus, we will work under the same general framework 
in order to better assess the applicability of our approach. In this sense, the usage of the 
classical version of the FDT is justified: 

where 

£(fc.c.) = - ''"^'''*-"' (3) 

UJ 

being the loss function of the system, which is assumed to be isotropic. The sum rules we 
employ are actually the power frequency moments of this loss function 

1 

C^[k) = - uj''C{k,uj)duj, z/ = 0,2,4. (4) 
Notice that odd-order moments vanish due to the symmetry o: 



" the loss function. 



The Nevanlinna formula of the classical theory of moments 



Ref. jlO|) expresses the response function 



(for a recent review see 



in terms of a Nevanlinna class function Q = Q{k, z), analytic in the upper half-plane Im 2; > 
and possessing there a positive imaginary part: lmQ{k,u + irj) > 0, t] > 0. The function 
Q{k,z) should also satisfy the limiting condition: 

^ 0, lmz>0. (6) 

Z zlcx) 



Any such function admits the integral representation 8|, |9| 

.00 / ^ t 



Q{k,z) = ih{k) + 



t- z 1 + ^2 



dg{t) 



(7) 



with Re h{k) > and some non-decreasing bounded function g (t) such that 

dg (t) 



< 00. 



The frequencies uJi{k) and uj2{k) in Eq. ([5]) are defined by the respective ratios of the 
moments Cy [U] 



the latter are expressible in terms of the system static characteristics, see Sect. Ill CI 

It is easily seen js] that the analytic prolongation of the loss function onto the upper 
half-plane Im z > constructed by means of the Cauchy integral formula, 



1 r 

C{k,z) = - 



ui — z 



admits the asymptotic expansion 
C{k, z — > 00) ~ 



Co (k) C2 (k) C4 (k) [ 1 

;: O 



(9) 



Z Z"^ z-^ \Z'-' 

while the expansion for the inverse dielectric function due to the Kramers-Kronig relations. 



e^{k,z) = l + -l ^-^—^du. 



(10) 



reads: 



e {k,z ^ 00) ^ 1 + 



C2{k) , C,{k) 



+ 



z^ 



Z'* 



1 + ^ + + o , 

z^ z^ \z^ 



Im 2; > 0. 



(11) 



Certainly, the expansions ([9]) and (fTTj) stem from ([5]) by virtue of the condition ([6]) [26|; the 
convergence of the moment Co (/c) follows from the existence of the static inverse dielectric 
function 

e-\k,^) = \mie-\k,%i]) = \^Py. H ^^^^^Llitl!^du;. (12) 

It is important that the inverse dielectric function ([5]) satisfies the sum rules ^ by 
construction, irrespectively of the form of the Nevanlinna parameter function Q[k,z) (see 



Sect. Ill Bp . In particular, this means that the asymptotic expansions ([9]) and (fTTl) are also 
valid for any function Q{k, z). 

Moreover, the inclusion of the moment Cn(k) guarantees that the inverse dielectric func- 
tion meets the compressibility sum rule 12] whenever e~^{k, 0) does. Within the present 
approach, this can be fulfilled by means of the relation between the charge-charge static 
structure factor and the static value of the plasma dielectric function (fT2|) . see below in 
Sect. ImXl Eqs. §^ and 

In general, there is no phenomenological basis for the choice of a unique Q{k,z), which 
would provide the exact expression for the loss function. The simplest is to approximate 
Q{k, z) by its static value ih (k) directly related to the static value Szz{k, 0) of the dynamic 
structure factor through Eq. ([2]), 

h(k) - ~ (131 

In this case the dynamic structure factor reads as 



S (kuj) = ^ - ^i') (u) 

' ^ 77/50 (k) ^2(^2 _ ^2)2 + /,2(^2 _ ^2^2 ' ^"^^^ 



This is the approximation which was used as a basis of the analysis carried out in Ref. 
Gj. If we use it now combined with more precise values for the power moments calculated 
using the Deutsch potential , we fail to predict the values of the Langmuir collective mode 
frequency, i.e., the position of the lateral peak of the dynamic structure factor, and its width, 
i.e., the approximation Q{k,z) = ih (k) is insufficient. 



B. The Nevanlinna parameter function 

It is quite clear from Eqs. ([2]), ([6]), and ([5]) that the parameter function Q{k, z) modifies 
both the position and the width of the Langmuir line in the spectrum of collective excitations 
of the system, as refiected in its charge-charge dynamic structure factor. 

If we neglect the processes of energy absorption completely, we should put in Eq. (IT^ 
h{k) = 0^. Mathematically this means the utilization, instead of the Nevanlinna formula 
([5]), of the canonical solution of the (truncated Hamburger) moment problem consisting of 
finding a loss function which satisfies the moment conditions @ or, equivalently, admits the 
expansion ([9]). 
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The canonical solution physically corresponds to the existence in the system of the dif- 
fusion (unshifted) and the infinitesimally-decaying Langmuir (shifted) mode, and implies a 
Feynman-like decay- free approximation for the charge-charge dynamic structure factor ([2]): 

,2 , ,2 , ,2 

]_ 

4 

To specify the Nevanlinna parameter function Q{k, z) and, hence, the non-canonical 
solution given by the Nevanlinna formula ([5]), we have to reconsider the details of energy 
absorption in the system without violating the sum rules. 

Precisely, in addition to the expansion (ITT]) , we want to satisfy the well-known Perel'- 
Eliashberg asymptotic form for the imaginary part of the dielectric function. The latter 
result can be summarized in the following way. 

In a completely ionized plasma for u ^ [I3h)'^ the microscopic acts of the electromag- 
netic field energy absorption become the processes which are inverse with respect to the 
bremsstrahlung during pair collisions of charged particles. 

As it was shown by L. Ginzburg [13], this circumstance permits to use the detailed 
equilibrium principle to express the imaginary part of the dielectric function, Ime {k^uo), of 
a completely ionized plasma in terms of the bremsstrahlung cross section and leads to the 
following asymptotic form of Im. e{k^uj) in a completely ionized (for simplicity, hydrogen- 
like) plasma obtained by Perel' and Eliashberg [14|] (see also Refs. [l^ and [ib] for an 
alternative derivation of this result and its specification based on the known expression for 
the bremsstrahlung differential cross section for high values of energy transfer and uj ^ 



Ime(fc,^»(/3n)-^)^^, (16) 
where 

Result (fT6!) also implies that higher even-order frequency moments, C21 (k) ,1 > 3, diverge. 
To take into account all convergent sum rules (power frequency moments) and the exact 
asymptotic relation (|T6l) we apply the Nevanlinna formula ([5]) with the following interpola- 



tion model expression for the Nevanlinna parameter function 271] 



Q{k,z) = B{k)^{l + t)+ th (k) = {k, z) + %Q2 (A;, z) , (18) 
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with 

AttA 

since by virtue of the Cauchy-Schwarz inequahty — ujf > 0. Observe that 



V2 r'J'u' 



p 



SO that 



35/4 ljI -LufV ujp^ ' ^TiVuf s{q, 0) ' 



UJ. 



n 



being the dimensionless charge-charge dynamic structure factor at = and q = ka, the 
values of s{q, 0) are provided in Table IV of Ref. [ij. 

The expression for the (inverse) dielectric function and, hence, for the dynamic structure 
factor ([2]) with the Nevanlinna parameter function determined in (fTSl) . leads to the correct 
static value Szz{k,0) of the dynamic structure factor, satisfies all three sum rules (jl]), and 
also satisfies the exact relation ( JT6l) . 

The range of frequencies studied in Ref. [l| was about (0,2c<jp), so that no data was 
obtained for the frequencies which satisfy the condition u ^ or, equivalently, 



UJ 

On the other hand, in spite of the classical approximations used here for comparison with Ref. 

as we have discussed previously, the system we consider possesses an inherent quantum 
mechanical nature, and we may presume that the asymptotic form ( fT6l) is applicable to it. 
By inserting Eq. f|T8|) into ([2]) we get 

g (k^^ _4 [ujl{k)-ujl{k)]Q2{k,u) 

^l3(l>ik)\u;[u;^-u;Uk)]+Q{k,u;)[u;^-u;f{k)f' ^ ' 

which reduces the determination of the dynamic structure factor to the knowledge of the 
static characteristics - the frequency moments Cy (k) , u = 0, 2, 4, see Sect. Ill CI 

The latter were calculated back in 1993 in Ref. Isj] in terms of the static structure 
factors of both system species. Sab {k), beyond the random-phase approximation with the 
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inclusion of both electronic and ionic local- field corrections (STLC). The electronic STLC 



was modelled 



18 



19| as an interpolation of the Geldart-Vosko form: 



(20) 



Cik'^ + C2k'^ ' 

where k^^ is the screening length of the electronic subsystem, treated as a one-component 
plasma (EOCP), an appropiate interpolation between the Debye and Thomas- Fermi radia 
(for details see Ref. [6|): 

ee (A;) = 1 + kl/e. (21) 



The interpolation parameters Ci an d Cg. were chosen in Ref. 19| to satisfy the compressibility 



sum rule for the classical EOCP 
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numerical simulations available at that time [20|, and the Kimball cusp-condition 



with the equation of state taken from the MD 

21: 



21 



The ionic STLC, 



lim Ge (k) 

fe— >oo 



G,{k) 



1 - ge (0) . 



(22) 



(23) 



cik"^ + C2k^ee (k) ' 

and the structural characteristics, i.e., the zero-separation value of the electronic radial 
distribution function, ge{0), and the static structure factors, Sab{k), were determined in Ref. 
6j according to the Ichimaru algorithm 23| in a self-consistent way in terms of the EOCP 
static structure factor Se{k). To take into account the quantum mechanical corrections, 
Se{k) was computed in Ref. using the temperature Green's function technique by a 
regularized summation over the Matsubara frequencies. 

To avoid all these cumbersome computations and, more important, to improve the agree- 
ment with the simulation data of Ref. [l| on the dynamic characteristics of dense two- 
component plasmas, the plasma static characteristics have been calculated in the present 
work by the HNC method with the Deutsch effective potential, as it was done in Ref. [l|. 



C. The moments 



The explicit form of the power moments C^(k) , u = 0,2,4, for the bare Coulomb potential 
I r I 

is known since long [llj] (for details see Ref. |l6l|): 



Co{k)=[l-e-'{k,0)] , C2{k)=ul, C^{k)=ul[l + W{k)] 



(24) 
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The moment Cq (k) , as it was already commented, is related to the static dielectric function 
of the system, the second moment is actually the /-sum rule, which is independent of the 
system interactions. The correction in the fourth moment contains different contributions: 

W{k) =K{k) + U{k) + H. 

The first contribution stems from the kinetic term of the system Hamiltonian. In the classical 
case, this coincides with the known Vlasov contribution to the dispersion relation: 

K{k) = 3^, (25) 

kj-, = Anne'^P being the square of the Debye wavenumber. Here, as in fl24|) . we only account 
for the electronic subsystem, due to the large asymmetry between the masses of electrons 
and ions. We use expression ( l25l) for comparison with the results of Ref. Nonetheless, 
due to the quantal nature of our system, it would be interesting to estimate quantitatively 
how the degeneracy would affect the dispersion law and the dynamic structure factor of the 
system through this kinetic contribution. In the quantum mechanical case, it can be recast 
as 

where the average of the square of the electron velocity is expressed as 

^ ' mp 

Fu{ri) = / Wl^^ 

Jo exp {x~r]) + l 

being the order- 1/ Fermi integral, and 1] = jSfi the dimensionless chemical potential of the 
electronic subsystem, which should be determined by the normalization condition 

The last two terms in the fourth moment correction term stem from the interaction 
contribution to the system Hamiltonian and are, therefore, dependent on the potential used. 
For the bare Coulomb potential we write: 

f oo 

27r% 



U{k) = -^l (See (P) - 1) / (P, k) dp, 







11 1 /"°° 

H = -K^ (0) = - ((7e,: (0) - 1) = / p'S,, ip) dp, 
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where we have introduced 



^ ~ 12 4P + 8pk^ 



In 



p + k 



p — k 



But for the sake of a better comparison with the MD results of Ref. [l|, we recalculated 
here these moments using the model Hamiltonian with the Coulomb potential substituted 
by the model Deutsch effective potential. Then, the moment C2 (k) (the /-sum rule) remains 
intact, the moment Co (k) changes together with the model system static dielectric function, 
and instead of C4 (fc) we have: 



C,{k) 



1 + W{k) 



(27) 



where the "model" W (k) has the same kinetic contribution K (k) (either classical or degen- 
erate), but the interaction contributions are substituted by 

1 



U{k) 



27r% 



P {See{p)-l) f (P, k) dp 



and 



Here 



H 



i^dp. 



p'^ + K 



(28) 



(29) 



2\2 



{P' + - k'f 
IQpk'^ 



In 



p + k 



In 



p — k 
{p + kf + 



{p - kf + < p2 + ^ 



.2 



In addition, the model partial static structure factors Sab (k) have been computed in 
the hypernetted-chain approximation using the Deutsch effective potential. This closes the 
algorithm of calculation of the static and dynamic characteristics of the system. 

Some numerical results are discussed in the following Section. 

III. NUMERICAL RESULTS 
A. Static characteristics 



As it was mentioned, we calculated the static structure factors and the radial distribution 
functions in the hypernetted-chains approximation using the Deutsch effective potential ([1]), 
just as it was done in Ref. l|. 
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Some data on the partial static structure factors are presented in Tables Hl lIIII It is not 



surprising that t 
Table El of Ref. 



le agreement we obtain with the values of the static characteristics given in 
ll, is within the computational error. We also add the corresponding values 



of the charge-charge static structure factor Szz{k) calculated by the following relation: 

SUk) = Su{k) + S,,{k) - 2SUk) (30) 

and check the zero-order frequency moment of the dynamic factor (fT9l) . 

1 

Szzik) = - j Szz{k,uj)duj; (31) 

to find an agreement to the fourth decimal digit. 

The other two dimensionless even-order power moments of the dynamic structure factor 
(fT9|) are defined as 

SAk) = / u'' Szzik, u;)du, z/ = 2,4, (32) 

the latter being provided in Table IIVI By virtue of the classical version of the FDT used 
here, the odd-order moments vanish due to the symmetry of f|T9|) . 

These values are then used to determine the characteristic frequencies ui (k) = 



6*2 (/c) /C^ik) and uj2{k) = y C4 (fc) /C2{k) which virtually coincide with their values 
calculated from the formulas (1271) . (l28l) . and (1291) . but now differ significantly from the 
values given in Table VI of Ref. [6] . 

Notice that due to (IT^ and the FDT, the static inverse dielectric function and the 
moment Co (k) are directly related to the charge-charge static structure factor (!3Tll : 

Szz (k) = -rrP.V- / Ime-i {k^uj) — = — (l - Ree-^ (k, 0)) , (33) 

J-00 '^^ I^D 

Thus the moment Cq {k) was estimated as 

Co{k) = ^Szz{k) (34) 

with the static structure factor Szz (k) also calculated in the hypernetted-chain approxima- 
tion using the Deutsch effective potential ([1]). 

Further, we display our results on the values of the partial radial distribution functions 
at zero separation, gee{0) and gie{0), computed using the effective potential ([1]) and also 
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taking into account (within the HNC scheme) the symmetry effects in the electron-electron 

exchange contribution to the effective potential (while leaving other components unchanged): 

, ,, ln2 / r^/tg \ 

^ee (r) = — (1 - exp (-Keer)) + — cxp ( "^^j^ I ; (35) 

for comparison we present also the values of 5'ie(0) calculated analytically in Ref. |2J|, see 
Table |Vl 

In addition, we display the graphs for the three partial radial distribution functions for 
the conditions F = 0.5 and = 0.4, calculated with the potential ([1]), Fig. [H The curves 
in Fig. [1] are virtually indistinguishable from those of Fig. 2 of Ref. [l| . 



B. Dynamic characteristics 

Results on the dynamic structure factor itself are broadly presented in Figs. [SHU 
The agreement between these graphs and the corresponding MD results of [1] is quan- 
titatively good, and it is better than that achieved in Ref. joj. The introduction of the 
non-constant Nevanlinna parameter function ffTSl) not only permits to obtain better agree- 
ment in the position of the Langmuir peaks, but also leads to the adequate broadening 
(damping) of the Langmuir mode. 

In this sense, it is interesting to calculate the complex solution for the dispersion equation 
e (/c, z) = explicitly in order to determine quantitatively the damping of the collective mode. 
From Eq. ffTTl) we get the equation 

z [z^ -uD+Q {z^ - ul) = 0. (36) 

Due to the fact that the function e~^{k,z) must be analytic in the upper half-plane, the 
solution of the dispersion equation, z{k), must possess a non-positive imaginary part, i.e., if 
z{k) = Rez(fc) + ilmz{k), then Im^; < 0. In particular, it is clear that for the Feynman-like 
approximation Q{h, z) = iO'^ we get for the (shifted) collective excitation the value of UJ2. 
However, for the model function f|T8l) . it is shown in Figs. [5] and |6] that Kez < 002, while the 
damping becomes more notorious. 

We observe that in all cases, except for F = 2 and the lowest values of g = ka, the 
effects of degeneracy produce stronger positive dispersion, i.e., the Langmuir mode frequency 
become higher than Up. In addition, degeneracy or quantum mechanical characteristics of 
the effective interaction produce also stronger damping of the mode. 
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Indeed, at higher frequencies, we deal with shorter distances or shorter times, where the 
non-Coulomb nature of the effective potential and the wave nature of electrons become more 
pronounced. 

Notice that these effects are reflected in the MD calculations using the Deutsch effective 
potential and are described in our calculations. 



IV. CONCLUSIONS 



The agreement with the MD data on the djTiamic structure factor and other dynamic 
characteristics of the model system, like the Langmuir collective mode dispersion, is improved 
with respect to the results obtained in Ref. [6|, and, simultaneously, important statistical 
characteristics of the model TCP in concern, like the probability to find both electron and 
ion at one point [g'ejlO)], are determined, the effects of direct and exchange interactions 
being compared. 

The importance of (l26l) and the applicability (ITSil are validated. 

Since the static characteristics of the system (the static structure factors and the radial 
distribution functions) were computed with high precision in the way employed in Ref. l|, 
the reliability of the calculation of the characteristic frequencies uji{k) and uj2{k) is improved. 
These quantities are essential to estimate the position and damping of the plasma collective 
mode. 

Another key ingredient for the good quantitative agreement achieved with the results of 
Ref. [l| with respect to the dispersion and decay of the Langmuir mode, is the introduction 
of a non-constant Nevanlinna parameter function accounting for the exact high-frequency 
asymptotic form of the system dielectric function. Still, further specification of this dynamic 
characteristic similar to the dynamic local-field correction in a TCP might be needed. 

We observed that for the Nevanlinna parameter function Q = Q{k, z) considered here, 
the value of the Langmuir frequency shifts from uJ2{k) closer to the plasma frequency, the 
effect of which might be considered correct from the experimental point of view. 

Further extension of the present mixed approach to more recent data for the model 
systems described by other effective potentials ^, 25| is planned. 
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Table Captions 



Table 1. Partial static structure factors at F = 0.5, = 0.4 with the Deutsch effective 
potential without exchange ([T]). 

Table II. Same as in Tabled but for F = 0.5, = 1. 

Table III. Same as in Tabled but for F = 2, = 1. 

Table IV. The fourth dimensionless power moment of Szzik^u), according to ([22D, 
with the Deutsch effective potential without exchange ([T]). In the classical case, the 
kinetic contribution is approximated by the Vlasov term fl25|) . whereas in the quantal case 
expression (126|) is used. 

Table V. Zero-separation values of the partial radial distribution functions compared to 



the results of Ref . [24 1 . 



18 



TABLE I: 



q — ka Sii[k) Si(.{k) S'ee(^) S^zik)"" 

0.767 0.5804 0.4387 0.6590 0.3620 

1.074 0.6257 0.3600 0.7391 0.6448 

1.381 0.6824 0.2813 0.8118 0.9316 

1.534 0.7118 0.2455 0.8425 1.0634 

"Calculated from ([30 l) or (|3T |) . 

TABLE II: 

q = ka Sii{k) Sie{k) See{k) Szz{kY 

0.767 0.6160 0.4606 0.6470 0.3418 

1.074 0.6663 0.3943 0.7144 0.5922 

1.381 0.7192 0.3275 0.7769 0.8412 

1.534 0.7447 0.2952 0.8081 0.9624 

''Calculated from (|30l) or ([3T |l . 
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TABLE III: 

q = ka Sii[k) Si(.{k) S'ee(/c) S^zik)"" 

0.767 0.5642 0.5821 0.7197 0.1198 

1.074 0.5133 0.5001 0.7385 0.2516 

1.381 0.5067 0.4275 0.7769 0.4286 

1.534 0.5174 0.3940 0.7993 0.5288 



"Calculated from ([30 l) or (|3T |) . 



TABLE IV: 







Classical 






Quantal 






r = 0.5 


r = 0.5 


r = 2.0 


r = 0.5 


r = 0.5 


r = 2.0 


= ka 


rs = 0.4 


rs = 1.0 


rs = 1.0 


rs = 0.4 


rs = 1.0 


rs = 1.0 


0.767 


0.8845 


0.9318 


0.1403 


1.1175 


0.9966 


0.1680 


1.074 


2.6028 


2.6943 


0.3294 


3.6028 


2.9856 


0.4462 


1.381 


6.2193 


6.3685 


0.6644 


9.3337 


7.3169 


1.0217 


1.534 


9.0728 


9.2555 


0.9073 


14.1572 


10.8366 


1.4855 
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TABLE V: 



r 


e 


b 






flpp(O)* 






0.1 


0.1000 


0.0184 


0.9598 


1.0684 


0.7309 


1.0474 


1.0358 


0.1 


2.0000 


0.3683 


0.7569 


1.5164 


0.3868 


1.5045 


1.3478 


0.5 


0.4344 


0.4000 


0.6629 


1.9865 


0.3897 


1.8874 


1.4412 


0.5 


1.0860 


1.0000 


0.4707 


3.5224 


0.2500 


3.2600 


1.8963 


1.0 


0.1000 


0.1842 


0.7738 


1.4979 


0.6159 


1.4276 


1.2351 



''With O. 

^With (Pi+dSSl). 

"^Values from Ref. 
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Figure Captions 



Fig. 1. The radial distribution functions Qabif), CL,b = e,i for F = 0.5, = 0.4, 
calculated with the potential ([T]). Solid lines are included to join the discrete points for a 
better visualization. 



Fig 2. Normalized djTiamic structure factor versus frequency as given by expression 
( IT9|) at F = 0.5, = 0.4. The model parameter function (fTSl) is used. The values of 
ka are: (1) 0.767, (2) 1.074, (3) 1.381, and (4) 1.534. (a) Classical calculations (solid 
lines) compared to the results of Ref. [l| with ka = 0.780 (diamonds), ka = 1.102 
(triangles), ka = 1.350 (boxes), and ka = 1.559 (stars), (b) Comparison between classical 
(dashed) and quantal calculations (solid). In the classical case, the kinetic contribution 
is approximated by the Vlasov term (125|) . whereas in the quantal case expression (126!) is used. 



Fig. 3. Same as in Fig. [2], but for F = 0.5, = 1. 



Fig. 4. Same as in Fig. [2|, but for F = 2, = 1. 



Fig. 5. Dispersion, z{k) = Ilez{k) + ilmz{k), for the collective excitation mode 
obtained from Eq. (|36l) in the classical approximation. The damped solution stemming 
from the model function (ITSj) is compared to the undamped solution corresponding to 
the Feynman-like approximation ( JT5l) . For expression (ITSil : F = 0.5, = 0.4 (triangles); 
F = 0.5, = 1 (stars); F = 2, = 1 (pentagons). For expression (fT5|) : F = 0.5, = 0.4 
(boxes); F = 0.5, = 1 (diamonds); F = 2, = 1 (circles). Sohd lines are included to join 
the discrete points for a better visualization. 



Fig. 6. Same as in Fig. [3, but for the quantal case. The degenerate kinetic term (1261) is 
used, instead of the classical one ( l25l) . 
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